Dual species matter qubit entangled with light 
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We propose and demonstrate an atomic qubit based on a cold 85 Rb- 87 Rb isotopic mixture, entan- 
gled with a frequency-encoded optical qubit. The interface of an atomic qubit with a single spatial 
light mode, and the ability to independently address the two atomic qubit states, should provide 
the basic element of an interferometrically robust quantum network. 
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Quantum mechanics permits the secure communica- 
tion of information between remote parties [l|, [2|, |3| . How- 
ever, direct optical fiber based quantum communication 
over distances greater than about 100 km is challenging 
due to intrinsic fiber losses. To overcome this limitation it 
is necessary to take advantage of quantum state storage 
at intermediate locations on the transmission channel. 
Interconversion of the information from light to matter 
to light is therefore essential. It was the necessity to 
interface photonic communication channels and storage 
elements that lead to the proposal of the quantum re- 
peater as an architecture for long-distance distribution 
of quantum information via qubits 0, [B| . 

Recently there has been rapid progress in interfacing 
photonic and stored atomic qubits. Two-ensemble encod- 
ing of matter qubits was used to achieve entanglement of 
photonic and atomic rubidium qubits and quantum state 
transfer from matter to light [6|. This was followed by 
a more robust single-ensemble qubit encoding [7J, which 
led to full light-matter-light qubit interconversion and 
entanglement of two remote atomic qubits Q. More 
recently, both two-ensemble and single-ensemble atomic 
qubits were reported using cesium gas 9, 10]. 

To realize scalable long distance qubit distribution 
telecommunication-wavelength photons and long-lived 
quantum memory elements are required [111 ]. Although 
multiplexing of atomic memory elements vastly improves 
the dependence of entanglement distribution on storage 
lifetime fl2| , there remains, the problem of robust atomic 
and photonic qubits for long-distance communication. 
Two-ensemble encoding suffers from the problem of long- 
term interferometric phase stability, while qubit states 
encoded in a single ensemble are hard to individually ad- 
dress. 

A protocol for implementing entanglement distribution 
with an atomic ensemble-based quantum repeater has 
been proposed It involves generating and transmit- 
ting each of the qubit basis states individually, in prac- 
tice via two interferometrically separate channels. Under 
prevailing conditions of low overall efficiencies it provides 
improved scaling compared to direct qubit entanglement 
distribution [J]. Its disadvantage is the necessity to sta- 
bilize the length of both transmission channels to a small 
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FIG. 1: (a) Schematic shows two orthogonal qubit states (ar- 
rows) encoded in two atomic ensembles coupled to distinct 
spatial light modes 0, @|, (b) the proposed architecture en- 
codes the states in a two species atomic mixture, coupled to 
a single light mode. 



fraction of the optical wavelength, as the distribution of 
qubit entanglement is sensitive to the relative phase fluc- 
tuations in the two arms. 

In this Letter we propose an interferometrically robust 
quantum repeater element based on entangled mixed 
species atomic, and frequency-encoded photonic, qubits, 
Fig. 1. This avoids the use of two interferometrically 
separate paths for qubit entanglement distribution. The 
qubit basis states are encoded as single spin wave excita- 
tions in each one of the two atomic species co-trapped in 
the same region of space. The spectroscopically resolved 
transitions enable individual addressing of the atomic 
species. Hence one may perform independent manipu- 
lations in the two repeater arms which share a single 
mode transmission channel. Phase stability is achieved 
by eliminating the relative ground state energy shifts of 
the co-trapped atomic species, as is in any case essential 
to successfully read out an atomic excitation [l3j |. 

We consider a co-trapped isotope mixture of 85 Rb and 
87 Rb, containing, respectively, N&5 and Ns7 atoms cooled 
in a magneto optical trap, as shown in Fig. 2. Unpo- 
larized atoms of isotope v (y S {85, 87}) are prepared 
in the ground hyperfine level where la*- 85 ') = 
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FIG. 2: Schematic of the experimental set-up showing the 
geometry of the addressing and scattered fields from the co- 
trapped isotope mixture of 85 Rb- 87 Rb. The write and read 
laser fields generate signal and idler fields, respectively de- 
tected at Dl and D2; Ei, E2 are optical frequency filters. 
PM1-4 are light phase modulators, <f) 3 and cj>i are relative 
phases of the driving rf fields, see text for details. 



I We consider the Raman configuration with 
ground levels \a^) and \b^) and excited level \c^) 



with energies Hui^ 



huj v \ and fiuj"' respectively. Level 



corresponds to the ground hyperfine level with 
smaller angular momentum, while level \c^ v 'j is the 



|5Pi/2) hyperfine level with Fc = F a l ■ A 150 ns long 
write laser pulse of wave vector k^, = k w y, horizontal 
polarization e# = z and temporal profile tp(t) (normal- 
ized to unity J dt \tp{t)\ 2 = 1) impinges on an electro- 
optic modulator (EOM), producing sidebands with fre- 
quencies cfcl 85 "* = ck w + 5u> w and cfcl 87 '' = ck w — 5u w 
(6uj w — 531.5 MHz) nearly resonant on the respective 
isotopic Di (\a^ v ') <-> c^)) transitions with detunings 

A„ = ck^ - {Jc ] - Ja ] ) w -10 MHz. Spontaneous 
Raman scattering of the write fields results in signal pho- 
tons with frequencies ck^ — ckw + — ui"' ) on the 
1 6^) <-+ | c ^) transitions. The positive frequency com- 
ponent of the detected signal electric field from isotope 
v with vertical polarization ey is given by 
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where u s (r ) is the transverse spatial profile of the 
signal field (normalized to unity in its transverse 
plane), and $ (t) is the annihilation operator for 
the signal field. These operators obey the usual free 
field, narrow bandwidth bosonic commutation relations 
[^" ) (t),^ )t (f)] = S vy S(t - f). The emission of V- 
polarizcd signal photons creates correlated atomic spin- 
wave excitations with annihilation operators given by 
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where 



cos 2 6 V = 
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X^/ a = C m ° m c C m b _ Q a ^ is a product of Clebsch- 
Gordan coefficients, and the spherical vector components 
of the spin wave are given by 
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The spin wave Zeeman components of isotope v are 
given in terms of the /i-th "Rb atom transition opera- 
tors u a {u) m - {v) m i and the write u w (r) and signal u s (r) 
field spatial profiles 
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The effective overlap of the write beam and the detected 
signal mode [3| is given by 
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where n^-'(r) is the number density of isotope za The 
interaction responsible for scattering into the collected 
signal mode is given by 

H s (t) = ih X (p(t) (cos ri^^it)^ 

+ smr 1 ^ 87 ^(t)s < - s7 ^) +h.c, (6) 



where \ = \f 
ramctcr, 



xi 5 



Xgy is a dimcnsionlcss interaction pa- 
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and crt^ are reduced matrix elements, riw is the 
average number of photons in the write pulse sideband 
with frequency ckw\ and the parametric mixing angle 
77 is given by cos 2 77 = xl 5 /(xl 5 + x| 7 )- The interac- 
tion picture Hamiltonian also includes terms represent- 
ing Rayleigh scattering and Raman scattering into unde- 
tected modes. One can show, however, that these terms 
commute with the signal Hamiltonian (Eq. ^) and with 
the operators ips (t) and s(") to order 0(1/ \fN). As a 
result, the interaction picture density operator for the 
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signal-spin wave system (tracing over undetected field 
modes) is given by UpoW , where po is the initial den- 
sity matrix of the unpolarizcd ensemble and the vacuum 
electromagnetic field, and the unitary operator U is given 
by 

\nU = xCcost?^ 85 )^ 85 ^ + sinr^ 87 )^ 87 ^ - h.c), (8) 

where d^ = J dtcp* (t)ips (t) is the discrete signal 
mode bosonic operator. When the write pulse is suffi- 
ciently weak we may write U — 1 = x(cos77d^ 85 ^s^ 85 ^ + 
sin?7d < - 87 ^,s < - 87 ^) + 0(x 2 ), i.e., the Raman scattering pro- 
duces entanglement between a two-mode field (frequency 
qubit) and the isotopic spin wave (dual species matter 
qubit). Although we explicitly treat isotopically distinct 
species, it is clear that the analysis is easily generalized 
to chemically distinct atoms and/or molecules. 

To characterize the nonclassical correlations of this sys- 
tem, the signal field is sent to an electro-optic phase mod- 
ulator (PM2 in Fig. 2) driven at a frequency Sui s — 
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)]/2 = 1368 MHz. 



The modulator combines the two signal frequency com- 
ponents into a central frequency ck s — c(/ci 85 ' ) + ki S7 ^)/2 
with a relative phase <fi s . A photoelectric detector pre- 
ceded by a filter (an optical cavity, El in Fig. 2) which 
reflects all but the central signal frequency is used to mea- 
sure the statistics of the signal. We describe the detected 
signal field using the bosonic field operator, 
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where e« S [0, 1] is the signal efficiency including prop- 
agation losses and losses to other frequency sidebands 
within PM2, and £i (t) represents concomitant vacuum 
noise. While quantum memory times in excess of 30 /is 
have been demonstrated (lBj . here the spin wave qubit 
is retrieved after 150 ns by shining a vertically polar- 
ized read pulse into a third electro-optic phase modulator 
(PM3 in Fig. 2), producing two sidebands with frequen- 

fo(85)) | c (85)) 



cies ckr and cki 87 ^ resonant on the 



and |?/ 87 )) <-> Ic^ 87 - 1 ) transitions, respectively. This re- 
sults in the transfer of the spin wave excitations to hor- 
izontally polarized idler photons emitted in the phase 



matched directions k = k w — k 



-k^. We treat 



the retrieval dynamics using the effective beam split 
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£r ir"\ where 



ter relations b^ u > = 
tr is the retrieval efficiency of the spin wave stored in 
the isotope "lib, U v) = J dt(p\ v) * (t)$ v) (t) is the dis- 
crete idler bosonic operator for an idler photon of fre- 
quency ckj" , <fi (t) is the temporal profile of an idler 
photon emitted from the y Rb spin wave (normalized 



to unity), and Tp^\t) is the annihilation operator for 
an idler photon emitted at time t. As with the signal 
operators, the idler field operators obey the usual free 
field, narrow bandwidth bosonic commutation relations 
[^ U) (t),ip^(t')} = 5 uy 5(t-t'). A fourth EOM, PM4, 
driven at a frequency Su>i — Su> w — (Ags + Ag7)/2 = 531.5 
MHz combines the idler frequency components into a 
sideband with frequency cki = c(fej 85 ' -I- fc l - 87 ' ) )/2 with 
a relative phase 4n . The combined idler field is measured 
by a photon counter preceded by a frequency filter (an 
optical cavity, E2 in Fig. 2) which only transmits fields 
of the central frequency cki. The detected idler field is 
described by the bosonic field operator, 



/ (85) / (87) 



+] j iz|L_ e */ 2 |f ) (t) + y iziL_ e -*/ 2 ef 

where ef^ G [0, 1] is the idler efficiency including prop- 
agation losses and losses to other frequency sidebands 
within PM4, and (i) represents associated vacuum 
noise. The write-read protocol in our experiment is re- 
peated 2 • 10 5 times per second. 

The signal-idler correlations result in phase- 
dependent coincidence rates given, up to de- 
tection efficiency factors, by C s j(0 s ,0j) = 

J dt s J dti($l(t s , <f>s)i>i(ti,<f)i)$ifc,<l>i)$a(ts,<t>a)\ 

From the state of the atom-signal system after the write 
process, UfoW, (EqJS}, we calculate the coincidence 
rates to second order in x, 

C si (cf> s , 4>i) = ^ (V 85) cos 2 r, + p^ sin 2 77 

+TvV 85 V 87) sin 277 cos (4>i - <t>s + 0o)) (9) 

where p^ = ef ej £a\ and T an d 0o represent a real 
amplitude and phase, respectively, such that 

fdt^ S5) *(t)^ 87 \t), (10) 
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and we account for classical phase noise in the rf driving 
of the EOM pairs PM1,4 and PM2,3, by treating (j) s and 
(pi as Gaussian random variables with variances 50 2 and 
<5</> 2 respectively, see Fig. 2. When the write fields are 
detuned such that the rates of correlated signal-idler coin- 
cidences are equal (i.e., when p^ 8 ^ cos 2 77 = p^- 87 ^ sin 2 77), 
the fringe visibility is maximized, and Eq. ([9]) reduces to 



= ^TM (85) cos 2 tj[1+T oos(&-^+0o)]. (11) 
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FIG. 3: Measured C s i{<j>a,4>i) as a function of ^ for <f) 3 — 
0, diamonds and for <f> 8 — —tt/2, circles. The angle <j>o is 
absorbed into the arbitrary definition of the origin, i.e., <j>o is 
defined to be zero. Solid lines are sinusoidal fringes based on 
Eq. JTTJ with T = 0.86. Single channel counts of Dl and D2 
show no dependence on the phases. 



Fig. 3 shows coincidence fringes as a function of (pi 
taken for two different values of <j> s . The detection rates 
measured separately for 85 Rb and 87 Rb were (a) 53 Hz 
and 62 Hz on Dl and (b) 95 Hz and 107 Hz on D2, respec- 
tively. These rates correspond to a level of random back- 
ground counts about 2.5 times lower than the minima of 
the interference fringes. This implies that the observed 
value of visibility T = 0.86 cannot be accounted for by 
random photoelectric coincidences alone. The additional 
reduction of visibility may be due to variations in the 
idler phases caused by temporal variations in the cloud 
densities during data accumulation, while the effects of 
rf phase noise are believed to be negligible. 

Following Ref. [l6| we calculate the correlation func- 
tion E(4> s ,(f>i), given by 



,^)-C 5i (^,^ i L )-C 5i (^,0 i ) + C 5 



C s 



,cf) l ) + C sl (<j> s ,(pl) + C sl (<t)^,(f) l ) + C. 



(12) 

where (pjr^ = (p s ^ + tt. We note that, by anal- 
ogy with polarization correlations, the detected signal 
[idler] field (t, (p^ ) is orthogonal to ^ s[i ] (t, (p s [i]), i.e., 

</>s[i]),i>l[i](t, <t>s[i]) = 0- 0ne nnds tnat a classi- 
cal local hidden variable theory yields the Bell inequality 
\S\ < 2, wtere S = ~ E(<P' s ,<Pi) ~ - 

E((/)' s , 0-) [13]. Using Eq. (fTT|) . the correlation function is 
given by 



E(<p s , fa) = T cos{(p s - (pi + (p ). 



(13) 



Choosing, e.g., the angles <p s — ~4>o, <Pi — 7r /4, <p' s — 
—(po — 7r/2, and <^ = 37r/4, we find the Bell parameter 
S = 2V2T. 

Table 1 presents measured values for the correlation 
function E((p s ,(pi) using the canonical set of angles 
cp s ,(pi- We find S exp = 2. 44 ±0.04 ^2 - a clear violation 



TABLE I: Measured correlation function E(<p a ,<pi) and S for 
At = 150 ns delay between write and read pulses; all the errors 
are based on the statistics of the photon counting events. 



i) 






-tt/2 

-7T/2 



tt/4 
3tt/4 
tt/4 
3tt/4 



0.629 ±0.018 
-0.591 ±0.018 
-0.614 ±0.018 
-0.608 ±0.018 
S exp = 2.44 ± 0.04 



of the Bell inequality. This value of S exp is consistent 
with the visibility of the fringes T w 0.86 shown in Fig. 
3. This agreement supports our observation that system- 
atic phase drifts are negligible. We emphasize that no 
active phase stabilization of any optical frequency field is 
employed. 

In conclusion, we report the first realization of a 
dual species matter qubit and its entanglement with a 
frequency-encoded photonic qubit. Although we em- 
ployed two different isotopes, our scheme should work for 
chemically different atoms (e.g., rubidium and cesium) 
and/or molecules. 
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